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Abstract. - We study the torque arising between two corrugated metallic plates due to the 
interaction with electromagnetic vacuum. This Casimir torque can be measured with torsion 
pendulum techniques for separation distances as large as 1/xm. It allows one to probe the 
nontrivial geometry dependence of the Casimir energy in a configuration which can be evaluated 
theoretically with accuracy. In the optimal experimental configuration, the commonly used 
proximity force approximation turns out to overestimate the torque by a factor 2 or larger. 



The Casimir effect [1] plays a major role in micro- and nano-electromechanical systems 
(MEMS and NEMS) [2]. Besides the normal Casimir force between metallic or dielectric 
plates [3], the observation of the lateral Casimir force between corrugated plates [4] opens 
novel possibilities of micro-mechanical control. The lateral force results from breaking the 
translational symmetry along directions parallel to the plates by imprinting periodic corru- 
gations on both metallic plates. As rotational symmetry is also broken by this geometry, a 
Casimir torque should also arise when the corrugations are not aligned. In the present letter, 
we study this effect which provides one with a new mechanism of micro-mechanical control to 
be exploited in the design of MEMS and NEMS. 

From the point of view of fundamental physics, this effect makes possible an accurate 
investigation of the non-trivial geometry dependence of the Casimir energy. Tests of this 
geometry dependence performed up to date have often been limited by the use of the proximity- 
force approximation (PFA) in the calculations [5]. Within this approximation, the energy is 
evaluated by simply averaging the Casimir energy between two planes over the distribution of 
local distances, which is certainly insufficient except for the limiting case of large corrugation 
wavelengths [6]. Other calculations [7] have gone farther than PFA while using a model of 
perfect reflection, which limited their validity to experiments performed with metallic mirrors 
at distances larger than a few /im. 

Here we use the scattering approach [8,9] extended to treat non-planar surfaces [10,11]. 
This perturbative approach is restricted by the assumption that corrugation amplitudes are 
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smaller than other length scales but it allows one to go beyond the range of validity of the 
PFA, with arbitrary values of the corrugation wavelengths with respect to plate separation 
and plasma wavelength of the metal. As shown below, it accurately describes configurations 
which are experimentally testable and where the Casimir energy has a non-trivial geometry 
dependence. This new effect induced by vacuum fluctuations may be small because the cor- 
rugation amplitudes are small, but the torque detection is particularly interesting from this 
point of view, since it allows one to use the exquisite sensitivity of torsion balances [12]. 

The idea of using torsion techniques has already been proposed to measure Casimir torques 
between anisotropic dielectric plates [13-16]. Recently, an experiment has been proposed to 
measure the torque with a birefringent disk on top of a barium titanate plate [17] . According to 
the authors, the proposed measurement is feasible at a plate separation L of lOOnm. We show 
below that the torque between corrugated metallic plates is up to three orders of magnitude 
larger than the torque between anisotropic dielectric plates, for comparable values of the 
separation distance and plate area. With realistic values for the corrugation amplitudes a\ 
and a 2 and wavelength Ac, we conclude that the torque measurement may be done at a 
distance L of the order of l//m. This minimizes the spurious effects due to the non-parallelism 
of the plates while making the contribution of interband transitions negligible. We thus model 
the finite conductivity of the metallic plates by the plasma model with a plasma wavelength 
Xp. We calculate the modification of the Casimir energy up to second order in the corrugation 
amplitudes assumed to be smaller than other length scales 



This condition, allowing us to use perturbation theory, will be satisfied by the numerical 
examples discussed below. 

The surface profile functions for the two plates are denoted as hj(r) with j = 1,2 labeling 
the two plates. They define the local heights, counted as positive when corresponding to 
separation decreases with respect to reference planes at z — and z — L. Both hj have zero 
spatial averages, so that L represents the average separation distance between the two surfaces, 
and r = (x, y) collects the two transverse coordinates defining the position on the plates. The 
corrugated surfaces are described by non-specular reflection coefficients that couple different 
field polarizations and momenta [18]. The second-order correction is then given by [11] 



where Hj{K) is the Fourier transform of hj(r). The response function Q(K) does not depend 
on the direction of the corrugation wavevector K. It is given in [11] as a function of the 
specular and non-specular reflection coefficients. Only the crossed terms of the form H1H2 
are kept in eq. ©, since terms quadratic in H± and in H2 do not depend on the relative angle 
between the plates, and hence do not contribute to the torque. 

We assume the corrugations to have sinusoidal shapes hj(r) = a.j cos(k 3 • r — kbj) with 
corrugation wavevectors kj having the same modulus k = 2tt/\c- The angle 9 between k x 
and k2 represents the angular mismatch between the two corrugations and k2 is supposed, 
for convenience, to be aligned along the direction of the x— axis. The parameters bj represent 
lateral displacements with respect to the configuration with a line of maximum height at the 
origin. With these conventions, eq. J5J) yields 



ai, a,2 •< L , Ac , Ap 



(1) 




(2) 



SEpp = ax Q{k) Re [, 



(3) 
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The integral yielding ^(ki) may be specified by considering that the corrugation h 2 is re- 
stricted to a rectangular section of area L x L y centered at x = 62, y = 0. We assume that both 
L x and L y are much larger than L, so that diffraction at the borders of the plates is negligible. 
The energy correction per unit area is then given by 

" L ' '' " " ! ° 2 Q(k) cos (kb) sinc{kL y sin0/2) ^ sine [kL x (l + e cos 9) /2] (4) 



L X Ly 2 



b = b 2 cos 9 — b\ , sinc(a;) = (5) 

x 

b is the relative lateral displacement along the direction of ki ; as expected by symmetry, the 
energy does not depend on displacements perpendicular to ki. Eq. Q satisfies reflection 
symmetry around the x and y directions (invariance under the transformations 9 — > —9 and 
9 — ► 7r — 9), which results from the fact that the corrugation lines have no orientation, ft 
contains as a special case (6 — 0) the result for pure lateral displacement derived in Ref. [If]. 

The energy variation with b and 9 given by eq. J1J is universal, since the separation distance 
L and additional parameters characterizing the metallic surfaces determine only the global 
pre- factor Q(k). In the limit of long corrugation lines, kL y 3> I, we may take the approximation 
9 <C f in eq. since 5Epp is negligible otherwise. If L x is smaller or of the order of L y , we 
find from eq. 

^pp = 0102 g ^ sinc ( kLy o/ 2 ). (6) 

L x L y 2 

Hence, the scale of variation of 9 is set by the parameter l/(kL y ) = Xc / (2tt L y ) . 

In Fig. ^ we plot SEpp (in arbitrary units) as a function of b and 9. Since Q(k) is neg- 
ative [11], the Casimir energy is minimum at 9 = and b = 0, Ac,2Ac, corresponding to 
the geometry with aligned corrugations and the local separation distance having the maxi- 
mum variation amplitude. There are also shallow wells around 9 ~ 1.43Ac / L y (minimum of 
sinc(fcLj / 0/2)) and b — Ac/2, 3Ac/2, ... If we start from 9 = b = and rotate plate 2 around 
its center, we follow the line b = in Fig. This figure shows that for small angles the plate 
is attracted back to 9 = b = without sliding laterally. On the other hand, if the plate is 
released after a rotation of 9 > Ac / L y its subsequent motion will be a combination of rotation 
and lateral displacement. In fact, the energy correction vanishes at 9 — Xc/L y , defining the 
range of stability of the configuration b = 9 = 0. Rotation is favored over lateral displacements 
for 9 < Xc/L y . 

We now proceed to an evaluation of the torque which is given by 

T=-^SE PP . (7) 
It is maximum at 9 = 0.66Ac / L y where it is given by 

— ^— = 0.109 ai a 2 kg{k)Ly. (8) 

L X Ly 

As could be expected, this maximal torque per unit area is proportional to the length L y of 
the corrugation lines, which provides the scale for the moment arm. In order to plot the value 
of the torque JSJ on Fig. we take the plasma wavelength corresponding to gold for both 
plates (Xp = 137nm). We chose corrugation amplitudes a\ — a 2 — 14 nm obeying condition 
and leading to a\a 2 — 200 nm 2 (to be compared with a\a 2 = 472 nm 2 in the lateral force 
experiment [4], where a\ and a 2 were unequal) and L y — 24 /im. Different values for a\a 2 and 
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Fig. 1 - Variation of Casimir energy (in arbitrary units) with the rotation angle 6 and the lateral 
displacement b. 

L y can also be derived from these results through a mere multiplication by suitable pre-factors 
(see eq. EJ, provided that these values satisfy ai,d2 <C L, Ac and L y ^> \c/{2ir). The first 
condition is at the heart of our perturbative approach, whereas the second condition can be 
relaxed by going back to the more general result (QJ. 

The dashed line in Fig. [21 corresponds to the corrugation period Ac = 1.2 /mi of the lateral 
force experiment [4]. At L = lOOnm, we find in this case t/(L x L v ) = 5.2 x lCP'N.m -1 , 
approximately three orders of magnitude larger than the torque per unit area for anisotropic 
plates calculated in Ref. [17] for the most favorable configuration at the same separation 
distance. The much larger figures found in our case should certainly allow one to perform 
the experiment at larger separation distances. For L ~ 1 /im, Fig. [3 shows that the torque 
for Ac = 1.2 /im starts to decrease exponentially. However, by selecting longer corrugation 
periods, one also finds measurable orders of magnitude in this range of distances. 

At any given value of L, the torque between corrugated plates can be made larger by 
chosing the corrugation period so as to maximize kQ(k). In the range of separation distances 
shown on Fig. |3 this corresponds to k ~ 2.6/L or Ac ~ 2irL/2.6. The torque is thus 
given by the dotted line, which provides upper bounds for its magnitude. At L = 1 /im, the 
optimum value is Ac = 2.4 /im, corresponding to the solid line in Fig. [21 In this case, we 
find t/(L x L v ) — 3.0 x 10~ 12 N.m _1 . As for the relevant angle scales, the stability threshold 
is at 9 = \c / L y = 0.1 rad = 5.7° and the maximum torque at 9 = 3.8°. If rotation is to be 
observed from the position of the tip of the rectangular plate, then the relevant parameter is 
the arc described by the tip (L y /2)9 ~ Ac, which is independent of L y . 

Under such optimum conditions, the torque probes a non-trivial geometry dependence of 
the Casimir energy, since the PFA thus grossly overestimates the effect, as shown now. The 
PFA holds for smooth surfaces corresponding to the limit Ac — > oo. It is recovered from 
eq. © at the limit k — > 0, where the response function satisfies the general condition [11] 
G(0) — ^pp(L), where epp is the Casimir energy per unit area for parallel planes. In this limit, 
the effect of geometry is trivial because the PFA directly connects the non-planar geometry to 
the more commonly studied configuration of parallel planes. We thus find a result determined 
by the Casimir energy for this simpler geometry: 




(9) 



In order to compare with the PFA, we plot our results for the torque as a function of k for 
L = 1 /on in Fig. (solid line), and with all other parameters as in Fig.|2 As discussed above, 
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Fig. 2 - Maximum torque per unit area as a function of mean separation distance. Corrugation 
amplitudes: aia2 = 200 nm 2 , line length: L v = 24 /im, plasma wavelength: \p — 137 nm. Solid line: 
Ac = 2.4 /im; dashed line: Ac = 1.2 /im; dotted line: Ac = 2irL/2.6 (optimum value). 



the torque is maximum at k = 2.6/ L — 2.6 /im . We also show the values obtained from the 
model with perfect reflectors (dashed line). They overestimate the torque by 16% near the 
peak region. 

According to eq. 10, the torque grows linearly with k in the PFA (dotted line). It is thus 
worth increasing the value of k and hence going out the region of validity of PFA. At the 
peak value k = 2.6/ L = 2.6 /im -1 , the PFA overestimates the torque by 103%. As the torque 
decays for larger values of k, the PFA rapidly becomes less and less accurate. The discrepancy 
is measured by the ratio 

m g( fc ) rim 
p(k) = = W) m 

At k = 0, the exact and PFA values coincide (p(0) = 1). For k > 0, p(k) < 1, as in 
the numerical example of Fig. 0] Thus, the PFA always overestimates the torque and the 
discrepancy increases with k, as expected, since smaller values of k correspond to smoother 
surfaces. 

For k ^> l/L, G{k) and p(k) decay exponentially to zero. Such behavior has a simple 




Fig. 3 - Maximum torque per unit area as a function of k — 2n/\c with L — 1 /im. Additional 
parameters are chosen as in Fig. 2. The solid line corresponds to the theory presented in this letter. 
We also plot the results for perfect reflectors (dashed line) and PFA with plasma model (dotted line). 
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Fig. 4 - Variation of p versus k. Parameters are chosen as in Fig. 3. Inset: variation of p(k) exp(fcL) 
(solid line) and perfectly-reflecting high-fc limit (dashed line). 

interpretation within the scattering approach [11]. The reflection by the corrugated surfaces 
produce lateral wavevector components of the order of k. After replacing the integral over real 
frequencies by an integral over the imaginary axis in the complex plane, we obtain propagation 
factors, representing one-way propagation between the plates, of the order of exp(— kL). 

In order to analyze the large-fc behavior in more detail, we plot p{k) exp(/cL) in the inset 
of Fig. 21 In the intermediate range 1/L <C k <C 2n/Xp = 46 fim , p(k) approaches the 
high-fc limit of the model with perfect reflectors [7], drawn as the dashed line in the inset of 
Fig.H 



As k approaches 2w/Xp, p(k) stays more and more below the result for perfect reflectors, 
reaching the following limit for very large values of k (not shown on the figure), 



Thus, the perfect-reflecting regime is reached only when the plasma wavelength is the shortest 
length scale in the problem (apart from the corrugation amplitudes). Just taking the limit 
L ^> Xp is not sufficient for recovering the limit of perfect reflectors [7]. This has again a 
very simple interpretation in the scattering approach. When L > Ap, the relevant input 
field modes are perfectly reflected because they have wavevectors of order of 1/L. However, 
diffraction by the corrugated surfaces produces wavevectors of the order of k, which are poorly 
reflected by the plates if k 3> 1/Ap. 

In conclusion, we have studied the torque induced by vacuum fluctuations between cor- 
rugated metallic plates. This Casimir torque may provide one with a new mechanism of 
micro- mechanical control to be exploited in the design of MEMS and NEMS. It could as well 
be a new observable of interest to test the geometry dependence of the Casimir energy We re- 
cover the PEA from our theory in the limit L <C Xq, but deviations become rapidly important 
as L/Xc is increased. 

The torque is up to three orders of magnitude larger than the torque between anisotropic 
dielectric plates for comparable distance and area. This should allow for an experimental 
observation of the Casimir torque at separation distances around L = 1 pm, using corrugation 



p(k) = — 7 (^-k) 4 exp(— kL) (perfect reflectors) 

77" 



(11) 




(12) 
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periods of the same order of magnitude (Ac = 2AL). In this configuration, the PFA grossly 
overestimates the torque by a factor of the order of 2. Moreover, it predicts an algebraic 
decay of the torque as L is increased past the optimum value, whereas the exact decay is 
actually exponential. This experiment would provide, for the first time, direct evidence of the 
nontrivial geometry dependence of the Casimir effect. 
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